Linearizer

ABSTRACT

The present invention is an improved linearizer that implements more complex transfer functions to provide the necessary linearization performance with a reasonable amount of signal processing resources. Particularly, the linearizer operates on an analog-to-digital converter and comprises a distortion compensator and one or more factored Volterra compensators, which may include a second-order factored Volterra compensator, a third-order factored Volterra compensator, and additional higher-order factored Volterra compensators. Inclusion of factored Volterra distortion compensators improves linearization processing performance while significantly reducing the computational complexity compared to a traditional Volterra-based compensator.

CROSS-REFERENCE TO RELATED APPLICATIONS

The present application is a continuation-in-part of and claims priority to co-pending U.S. patent application Ser. No. 13/455,071, filed Apr. 24, 2012, and entitled “Adaptive Digital Receiver”; which is a continuation-in-part of U.S. patent application Ser. No. 12/908,783, filed Oct. 10, 2010, and entitled “Mismatch Compensators and Methods for Mismatch Compensation” (now U.S. Pat. No. 8,164,496); which is a continuation of U.S. patent application Ser. No. 12/817,075, filed Jun. 16, 2010, and entitled “Amplifier Linearizer” (now U.S. Pat. No. 7,940,198); which is a continuation-in-part of U.S. patent application Ser. No. 12/112,380, filed Apr. 30, 2008, and entitled “Adaptive Mismatch Compensators and Methods for Mismatch Compensation” (now U.S. Pat. No. 7,782,235), the disclosures of which are all herein incorporated by reference in their entirety. The present application also claims priority to U.S. Provisional Patent Application No. 61/658,116, filed on Jun. 11, 2012, and entitled “Linearity Compensator”; the disclosure of which is herein incorporated by reference in its entirety.

BACKGROUND OF THE INVENTION

1. Field of Invention

The present invention relates to electronics and more specifically, to techniques for improving the linear performance of electrical components such as analog-to-digital converters and power amplifiers in radio frequency (RF) transceiver systems.

2. Description of Related Art

An electronic or power amplifier is a device for increasing the power of an input signal. Power amplifiers in most electronic systems are required to be “linear,” in that they must accurately reproduce the signal present at their input to achieve efficiency and spectral purity. An amplifier that compresses its input or has a non-linear input/output relationship causes the output signal to splatter onto adjacent frequencies. This causes interference on other frequency channels. Power amplifiers tend to become more non-linear as their power increases towards their maximum rated output.

Power amplifier linearity and efficiency are crucial issues in the design of electronic systems. The power amplifier is one of the most import parts of and usually the largest single contributor to the power consumption of an RF system. Thus, it is desirable to maximize the efficiency of a power amplifier. However, the more efficient power amplifier configuration is used the more nonlinear it usually becomes. Linearity errors in power amplifiers cause harmonic distortion and intermodulation distortion, which can limit the performance of state-of-the-art electronic systems such as, but not limited to radar systems, digital transceivers for wireless 3G and 4 G communications, laboratory test equipment, medical imaging, and audio and video compression. Reducing errors in digital-to-analog converters, analog-to-digital converters, sample-and-hold circuitry, and buffer and power amplifiers can significantly improve the performance of the critical conversion process.

An analog-to-digital converter (ADC) is a device for converting continuous radio frequency signals into discrete-time sampled, quantized data for subsequent digital processing. Like power amplifiers, analog-to-digital converters in most electronic systems are required to be “linear,” in that they must accurately reproduce the signal present at their input to provide a high-resolution digitized version at its output. An analog-to-digital converter with a non-linear transfer function will introduce distortion components, such as harmonic or intermodulation distortion, that limit the effective resolution and dynamic range of the analog-to-digital conversion.

Several common linearization methods exist for improving the linearity of devices such as power amplifiers and analog-to-digital converters, all of which suffer from drawbacks. All conventional linearization methods are limited in their maximum correctable range, which is the region of power output level near the onset of saturation. One method, known as a feedforward technique, is frequently employed and improves linearity, but results in poor power amplifier efficiency.

Predistortion is another technique used to improve the linearity of amplifiers. Digital predistortion is used to linearize the nonlinear response of a power amplifier over its intended power range. A predistortion circuit inversely models an amplifier's gain and phase characteristics and, when combined with the amplifier, produces an overall system that is more linear and reduces the amplifier's distortion. In essence, “inverse distortion” is introduced into the input of the amplifier, thereby cancelling to some degree any non-linearity the amplifier might have. However, the effectiveness of any predistortion technique is directly dependent on the accuracy of the predistortion transfer function, i.e., the model of the amplifier's gain and phase distortion characteristics. Traditional predistortion techniques implement a second-order or third-order polynomial as the transfer function, which improves the linearity of a power amplifier. However, for advanced RF systems with very high instantaneous bandwidths, transfer functions based on second-order or third-order polynomials are not accurate enough to prevent all non-linear distortion. In fact, many RF devices produce irregular nonlinearities that are difficult to model with standard, integer-power polynomial functions. Moreover, traditional predistortion techniques typically employ only one transfer function for a power amplifier, which may be suitable for a particular set of operating conditions. However, when operating conditions vary, e.g., temperature, time, or frequency, the transfer function may no longer be suitable to adequately remove non-linear distortion. Accordingly, there is a need for an improved linearization technique that implements higher-order transfer functions and accounts for varying operating conditions.

Digital post-processing is another technique used to improve the linearity of analog-to-digital converters. Digital post-processing is used to linearize the nonlinear response. Digital post-processing inversely models an analog-to-digital converter's nonlinear distortion transfer function such that the deleterious nonlinear distortion components can be subtracted from the output of the analog-to-digital converter to produce an overall system that is more linear and reduces the nonlinear distortion. In essence, “inverse distortion” is introduced into the output of the analog-to-digital converter, thereby cancelling to some degree any non-linearity the analog-to-digital converter might have. However, the effectiveness of any digital post-processing technique is directly dependent on the accuracy of the nonlinear distortion transfer function, i.e., the model of the analog-to-digital converter's nonlinear distortion characteristics. Traditional digital post-processing techniques implement a second-order or third-order polynomial as the transfer function, which improves the linearity of an analog-to-digital converter. However, for advanced RF systems with very high instantaneous bandwidths, transfer functions based on second-order or third-order polynomials are not accurate enough to prevent all non-linear distortion. Traditional higher-order or more complex models, such as Volterra-based compensation techniques, require a prohibitive amount of signal processing resources (such as multiply-accumulate functions which increases the physical size, weight, power, and cost of the hardware implementation of the digital post-processing). Accordingly, there is a need for an improved linearization technique that implements more complex transfer functions to provide the necessary linearization performance with a reasonable amount of signal processing resources.

SUMMARY OF THE INVENTION

The present invention overcomes these and other deficiencies of the prior art by providing a linearizer comprising a distortion compensator and one or more factored Volterra compensators, which may include a second-order factored Volterra compensator, a third-order factored Volterra compensator, and additional higher-order factored Volterra compensators. Inclusion of factored Volterra distortion compensators improves linearization processing performance while significantly reducing the computational complexity compared to a traditional Volterra-based compensator.

In an embodiment of the invention, a linearity compensator for removing nonlinear distortion introduced by an electronic system comprises: a distortion compensator; and one or more factored Volterra compensators. The one or more factored Volterra compensators comprise a second-order factored Volterra compensator, a third-order factored Volterra compensator, and/or additional higher-order factored Volterra compensators, wherein the additional higher-order factored Volterra compensators are fourth or higher-order factored Volterra compensators. The second-order factored Volterra compensator comprises a number of parallel finite impulse response filters and a number of memoryless second order polynomials. The second-order factored Volterra compensator is implemented in a parallel polyphase configuration. The third-order factored Volterra compensator comprises a number of parallel finite impulse response filters and a number of memoryless third order polynomials. The third-order factored Volterra compensator is implemented in a parallel polyphase configuration. The electronic system comprises an analog-to-digital converter or a power amplifier. The Volterra compensators comprise a number of linear filters implementing factored Volterra kernels. The factored Volterra kernels are measured via harmonic probing comprising multi-tone test signals injected into the input of the electronic system. The number tones in the multi-tone test signals is greater than or equal to the order of the corresponding factored Volterra kernel. Alternatively, the factored Volterra kernels are measured via harmonic probing comprising random test tone signals with non-overlapping distortion components. The linearity compensator comprises an adaptive nonlinear distortion estimator. The factored Volterra kernels are calculated using an over-constrained linear least mean square filter design. Alternatively, the factored Volterra kernels are calculated using interpolation and extrapolation and inverse fast Fourier transforms. The distortion compensator comprises: a first linear filter providing a fundamental signal path and a compensator signal path; a distortion transfer function operating on the compensator signal path, wherein the distortion transfer function comprises a nonlinear polynomial and a second linear filter; and an adder for summing the fundamental signal path and the compensator signal path. The linearity compensator may further comprise a bandpass interpolator operating on the compensator signal path. The bandpass interpolator comprises an upsampler and a bandpass filter. The linearity compensator may further comprise a downsampler operating on the compensator signal path. The filter coefficients of the first linear filter are measured by fitting the filter coefficients to a coherently measured amplitude and phase of a transfer function of the distortion transfer function. The distortion transfer function further comprises an adaptive nonlinear distortion estimator. The filter coefficients of the first linear filter are iteratively optimized. The filter coefficients of the first linear filter are optimized using MINIMAX optimization. The filter coefficients of the second linear filter are optimized using MINIMAX optimization. The factored Volterra compensators may comprise a subset of all the factors.

The present invention provides numerous advantages such as improving and maintaining extremely high performance of power amplifiers and analog-to-digital converters over a wider range of signal and environmental conditions, with reasonable amount of signal processing resources and adaptively tracking and correcting parameters which may drift with time, temperature, frequency, amplitude, etc. The present invention simplifies or eliminates outboard calibration, significantly improves power amplifier and analog-to-digital converter performance, provides built-in test capability (i.e., system failures can be detected and analyzed), and increases reliability. The adaptive processing algorithms of the present invention efficiently share processing resources already present in adaptive systems (e.g., adaptive digital beamforming), and reduce the performance requirements of the analog circuitry because the digital signal processing is used to dramatically improve performance. In sum, the present invention lowers the cost, power, and size of power amplifiers and analog-to-digital converters while improving linearity and efficiency.

The foregoing, and other features and advantages of the invention, will be apparent from the following, more particular description of the preferred embodiments of the invention, the accompanying drawings, and the claims.

BRIEF DESCRIPTION OF THE DRAWINGS

For a more complete understanding of the present invention, the objects and advantages thereof, reference is now made to the ensuing descriptions taken in connection with the accompanying drawings briefly described as follows:

FIG. 1 illustrates a block diagram of a distortion compensator;

FIG. 2 illustrates a block diagram of the distortion transfer function circuitry of FIG. 1;

FIG. 3 illustrates exemplary asymmetric and symmetric distortion transfer functions;

FIG. 4 illustrates an adaptive linearized power amplifier system;

FIG. 5 illustrates the adaptive nonlinear distortion estimator of FIG. 1;

FIG. 6 illustrates a heuristic calibration system;

FIG. 7 illustrates a process for heuristically compensating nonlinear distortion;

FIG. 8 illustrates a block diagram of a linearity compensator according to an embodiment of the invention;

FIG. 9 illustrates the structure of the second-order factored Volterra compensator according to an embodiment of the invention;

FIG. 10 illustrates the second-order factored Volterra compensator operating in a parallel polyphase configuration according to an exemplary embodiment of the invention;

FIG. 11 illustrates filters implemented in a parallel polyphase configuration according to an embodiment of the invention;

FIG. 12 illustrates the structure of the third-order factored Volterra compensator according to an embodiment of the invention; and

FIG. 13 depicts a distortion compensator where a subset of the full Volterra kernel is implemented in a Linear-Nonlinear-Linear (LNL) Hammerstein-Wiener non-linear error compensation model.

DETAILED DESCRIPTION OF EMBODIMENTS

Further features and advantages of the invention, as well as the structure and operation of various embodiments of the invention, are described in detail below with reference to the accompanying FIGS. 1-13, wherein like reference numerals refer to like elements. Although the invention is described in the context of power amplifiers and analog-to-digital converters, one of ordinary skill in the art readily appreciates that the linearization techniques described herein are applicable to any type of electronic component where it is desired to accurately and adequately eliminate inherent non-linear errors.

In general, predistortion linearization circuitry implements, among other things, (1) a higher-order polynomial model of an amplifier's gain and phase characteristics—higher than a third-order polynomial model; (2) an adaptive calibration technique; and (3) a heuristic calibration technique. The higher-order polynomial model is generated by introducing, for example, a plurality of multi-tone test signals with varying center frequency and spacing into the power amplifier. From the power amplifier's corresponding output, the nonlinearities are modeled by employing a higher-order curve fit to capture the irregularities in the nonlinear transfer function. Different transfer functions can be implemented for different operating conditions. Adaptive calibration is based on a feedback analysis technique, which updates the applicable transfer function by analyzing the error signal between the introduced uncompensated input signal and the compensated output signal in real-time. Heuristic calibration implements different transfer functions based on historical operating conditions and optimal configurations of the power amplifier.

FIG. 1 illustrates a block diagram of a distortion compensator 100. The distortion compensator 100 operates in the digital domain and serves as a general linearity error compensator that accurately nulls distortion by combining a predistortion signal e[n] with an input signal y[n] to produce a compensated output signal z[n], such that:

y[n]+e[n]=z[n]

The distortion compensator 100 is coupled to an RF transmission chain (not shown) and digital signal processing circuitry (not shown) as part of an RF communications system. A typical RF chain comprises a digital-to-analog converter, one or more filters, upconversion circuitry, a high-power amplifier (HPA), and an RF antenna. Although the present disclosure is described in the context of predistortion, i.e., pre-processing circuitry, the linearization techniques described herein can be implemented as digital post-processing circuitry to compensate for distortion introduced by an RF receiving chain. A typical RF receiving chain (not shown) comprises an RF antenna, a low-noise amplifier (LNA), downconversion circuitry, one or more filters, and an analog-to-digital converter. The implementation of an RF transmission chain, RF receiving chain, and digital signal processing circuitry in an RF transceiver system is readily apparent to one of ordinary skill in the art.

The distortion compensator 100 comprises an upsampler 110, a band-pass filter 120, a distortion transfer function circuitry 130, a downsampler 140, and an adaptive nonlinear distortion estimator 150. The inclusion of the upsampler 110, the band-pass filter 120, and the downsampler 140 is optionally employed in order to perform interpolation, the implementation of which is apparent to one of ordinary skill in the art. Interpolation increases the bandwidth of the distortion compensator 100 by estimating the compensator input signal, y[n], at intermediate points. This allows the compensator 100 to properly resolve harmonic and intermodulation distortion components that exceed the Nyquist bandwidth of the sampled system. In operation, the upsampler 110 increases the sampling rate of the compensator input signal, y[n], by a factor of M. The band-pass filter 120 preserves the fundamental content of the compensator input signal, y[n], by removing content above and/or below the original Nyquist limit of the signal. The downsampler 140 decreases the sampling rate to the original sampling rate of the compensator input, y[n]. The downsampler 140 is coupled to an adder as shown. Optionally, a linear filter may be included on the fundamental path which operates on the compensator input signal, y[n], and the output of which couples to the adder instead of directly coupling the compensator input signal, y[n], to the adder as shown. This optional linear filter can be configured to compensate for the linear, frequency-dependent gain and phase response, which when compensated, also improves the frequency-dependent performance of the nonlinear distortion.

FIG. 2 illustrates a block diagram of the distortion transfer function circuitry 130. The distortion transfer function circuitry 130 is coupled to the band-pass filter 120 and the downsampler 140, and receives input signals 155 and 165 from the adaptive nonlinear distortion estimator 150. The distortion transfer function circuitry 130 can employ any number “N” of nonlinear distortion models. Referring to FIG. 2, the distortion function circuitry 130 comprises N number of nonlinear distortion models 210A-N and an equal number of gain/phase shift filters 220A-N. Each nonlinear distortion model 210I is coupled to a corresponding gain/phase shift filter 220I. Each gain/phase shift filter 220 comprises one or more digital filters (e.g., finite impulse response (FIR) filters of length 16 taps). In an embodiment of the invention, the digital filters are fit to the measured amplitude and phase variation across the desired frequency band using standard least mean squares (LMS) or minimization of maximum error (MINIMAX) filter design algorithms. The digital filters are key components that allow precise cancellation of the nonlinear distortion components as they vary in gain and phase over frequency.

Two nonlinear distortion models 210A and 210B can be implemented. The nonlinear model 210A compensates for asymmetric, or odd order, distortion components. The nonlinear model 210B compensates for symmetric, or even order, distortion components. The gain/phase shift filter 220A or 220B is coupled to the respective nonlinear model 210A or 220B to dynamically adjust the distortion components such that they are canceled in the compensated output z[n] over a wide frequency band as described below.

Each nonlinear distortion model 210 is produced by curve fitting a measured distortion transfer function (i.e., distortion level as a function of input level). For example, curve fitting may be implemented using a MINIMAX technique, the implementation of which is apparent to one of ordinary skill in the art. Other curve fitting techniques can be employed such as, but not limited to an LMS methodology or a nonlinear statistical regression analysis, the implementation of all of which are apparent to one of ordinary skill in the art. In at least one embodiment of the invention, the nonlinear distortion models 210A and 210B comprise a 13^(th) order polynomial as the asymmetric transfer function and a 12^(th) order polynomial as the symmetric transfer function. However, higher order polynomials can be implemented if adequate processing resources are available. Nonetheless, from experimental results, 7^(th) order and 8^(th) order polynomials are usually adequate to compensate for asymmetric and symmetric non-linear distortion when practicing the present invention. In another embodiment of the invention, a single nonlinear distortion model 210 is employed by implementing an even and odd order polynomial model. Alternatively, a single polynomial model is implemented where one or more of the exponents are not integers.

The asymmetric and symmetric models 210A and 210B are fitted to an error distortion profile resulting from a plurality of injected test signals. For example, a set of 10 to 15 multi-tone test signals with varying center frequency and spacing is used to measure the resulting harmonic and intermodulation distortion components likely during manufacturing of the RF system. For each test signal, a buffer of approximately 128K samples is captured and the fast Fourier transformation (FFT) spectrum is analyzed, the implementation of which is apparent to one of ordinary skill in the art, to measure the relative amplitude and phase shift of each of the distortion components. The results are tabulated, and distortion transfer functions are calculated for the selected distortion components. The distortion transfer function is the measured level of the selected distortion components as a function of the input signal level. A higher order polynomial function as described above is fit to the measured distortion transfer functions to form the nonlinear distortion model, which represents the typical distortion behavior of the device.

FIG. 3 illustrates exemplary asymmetric and symmetric distortion transfer functions 300 and 310. The asymmetric distortion function 300 accounts for odd order distortion components (e.g., up to the seventh-order) that are not possible to model accurately with a standard third-order polynomial function 305. The symmetric distortion transfer function 310 accounts for even order distortion components (e.g., up to the eighth-order) that are not possible to model accurately with a standard second-order polynomial function 315.

Referring back to FIG. 2, in an alternate embodiment of the invention, the nonlinear models 210A-N are configured for non-harmonically related distortion components, for example, mixing image products, clock leakage, adjacent channel leakage, and/or other interference or nonlinear distortion. To accurately model mixing image products, the nonlinear model 210A-N may include a mixer, which multiplies a signal with a sinusoidal local oscillator signal to shift its frequency content to another frequency. It may optionally include a bandpass filter to select one or more mixing image products.

FIG. 4 illustrates an adaptive linearized power amplifier system 400. The system 400 comprises the compensator 100, an RF transmitter 420, a power amplifier 410, and an RF receiver 430. The RF transmitter 420 comprises a digital-to-analog converter and may optionally include RF upconversion electronics to convert a baseband signal to a higher RF. The power amplifier 410 amplifies the RF signal and introduces the undesired nonlinear distortion. The power amplifier 410 output is tapped and coupled to the RF receiver 430. The RF receiver 430 comprises an analog-to-digital converter and may optionally include RF downconversion electronics to convert the power amplifier output from RF to baseband. The feedback signal 465 is a digitized representation of the tapped power amplifier output, which is used by the compensator 100 to update its behavior, through the adaptive nonlinear distortion estimator 150, for optimal performance.

FIG. 5 illustrates the adaptive nonlinear distortion estimator 150. The estimator 150 comprises a delay 510 and a parameter updater 520. The delay 510 delays the system input, y[n], so that it properly aligns with the feedback signal 465, which is a digitized representation of the power amplifier 410 output and is subject to delays incurred in the processing and circuitry of the system 400. The delayed input signal is combined with the feedback signal 465 via an adder to subtract out an error signal between the compensated output signal, z[n], and the uncompensated signal, y[n]. The error signal is a representation of the difference between the input uncompensated signal, y[n], and the amplifier output, which includes the nonlinear distortion components that the compensator 100 corrects.

This error signal is introduced to the parameter updater 520. The parameter updater 520 comprises a spectrum analyzer 530, a channelizer 540, a transfer function calculator 550, a coefficient generator 560, and a controller 570, the implementation of all of which are apparent to one of ordinary skill in the art. Upon processing of the error signal, the transfer function calculator 550 outputs updated transfer function parameters 155, which are used by the nonlinear distortion transfer function circuitry 130 to further refine the nonlinear models 210A-N. The coefficient generator 560 outputs updated filter coefficients 165, which are used by the gain/phase shift filters 220A-N to fine tune the overall predistortion. The spectrum analyzer 530 evaluates the spectral content of the error signal, for example, with an FFT operation. The controller 570 uses this spectral content information to determine if the current signal is suitable for estimation of the nonlinear distortion components. For example, if the signal levels are currently too low for accurate estimation, then the controller will wait for a larger signal to appear before continuing with the estimation. The channelizer 540 partitions the error signal into frequency subbands, for example, with a digital filter bank, to analyze the distortion components over different frequencies.

The transfer function calculator 550 stores a table of the current distortion transfer function that maps the delayed input signal to the error signal, for example, in random access memory (RAM). The transfer function calculator 550 may also use a curve-fitting algorithm, the identification and implementation of which are apparent to one of ordinary skill in the art, to fit a polynomial equation to the transfer function. The calculator 550 outputs transfer function parameters 155, which may be this equation or the memory table. The calculator 550 optionally outputs similar transfer function parameters 555, which may be this equation or the memory table, for use in a heuristic compensation process described below. The coefficient generator 560 uses the output of the spectrum analyzer 530 to measure the relative gain and phase shift of the selected nonlinear distortion components. These gain and phase shift measurements are cataloged over frequency for different input signals. Standard digital filter design methods may be used to fit a digital filter to the gain and phase shift measurements to generate the coefficients 165 of the gain/phase shift filters 220A-N. The generator 560 optionally outputs coefficients 565 for use in a heuristic compensation process described below.

In essence, the adaptive nonlinear distortion estimator 150 serves as a blind adaptive calibration routine to significantly simplify or even eliminate system calibration, i.e., the adaptive routine can be used to monitor the system output, z[n], and interactively adjust the nonlinear models 210A-N and the gain/phase shift filters 220A-N to minimize the nonlinear distortion without requiring any external calibration signal sources or by requiring the interruption of the normal operation of the system. The blind adaptive calibration routine does not require the use of a pilot signal to characterize the system, which eliminates the need to generate test signals during operation of the applicable RF communications system. Since the mathematical relationship between the desired output and the nonlinear error signals is known for arbitrary input signals, this relationship can be used to estimate the amplitude and phase of the error signals. The blind adaptive calibration routine (implemented, for example, in a digital signal process) iteratively adjusts the digital filters to subtract (or null) these error signals from the output signal. FFT analysis (or other wavelet time-frequency analysis) can be used to measure how accurately the error signals have been nulled. Accurate nulling of these error signals means the processing is optimized.

Heuristic calibration is also implemented. Heuristic calibration maintains and updates the calibration history over numerous operating conditions, e.g., time, temperature, frequency, etc. For example, distortion transfer functions are generated for various sets of operating conditions and are stored in a look-up table. In the event that operating conditions change, an applicable higher order asymmetric and symmetric transfer function can be loaded into the nonlinear models 210A-N or an applicable set of filter coefficients 165 can be loaded into the gain/phase shift filters 220A-N, or both. This permits the compensator 100 to quickly switch to an optimal configuration, which will be further refined through the blind adaptive calibration technique discussed above.

FIG. 6 illustrates a heuristic calibration system 600. Particularly, the heuristic calibration system 600 comprises a nonlinear distortion cataloger 610, which stores a history of various measurements, including, but not limited to transfer function parameters 555, filter coefficients 565, and other relevant parameters such as mode of operation 635, temperature 645, and time 655. The nonlinear distortion cataloger 610 outputs distortion transfer function parameters 155 and updated filter coefficients 165 that are calculated heuristically based on the current operation conditions, such as time, temperature, and mode of operation. The nonlinear distortion cataloger 610 can record a number of different measurements to allow for accurate interpolation or extrapolation of the transfer function parameters 155 and updated filter coefficients 165, given the RF system's current operating mode, current temperature, other current conditions, and over a range of frequencies. This record of different measurements can be stored in a catalog (not shown). The nonlinear distortion cataloger 610 can extract the most relevant measurements from the catalog in response to the current operating conditions of the system.

For example, the nonlinear distortion cataloger 61 can extract all of the nonlinear distortion level measurements that have been stored in the catalog that are for the current operating mode 635, within 5 degrees Celsius of the current temperature 645, with the last 30 seconds of the current time 655, and over all the frequencies at which measurements have been stored in the mismatch catalog. Repeated measurements at the same frequency can be averaged for a more accurate estimation. Missing measurements for particular frequencies can be interpolated or extrapolated from the extracted measurements. Measurements that vary significantly from previously stored measurements can be discarded since they may indicate a measurement glitch. This process performed by the nonlinear distortion cataloger 610 of extracting the relevant data from the catalog, combining repeated measurements, interpolating or extrapolating missing measurements, and discarding inaccurate measurements is a heuristic process that generates the transfer function parameters 155 and updated filter coefficients 165.

FIG. 7 illustrates a process 700 for heuristically compensating nonlinear distortion. This process 700 is controlled by the nonlinear distortion cataloger 610 as shown in FIG. 6. The nonlinear distortion cataloger 610 can be implemented as a microprocessor or an embedded processor with a field programmable gate array. Particularly, the process 700 begins by evaluating (step 705) the amplitudes and frequency content of the nonlinear distortion levels and determining (step 710) if that data is suitable for accurate estimation of the nonlinear distortion. The accuracy of the estimation can be degraded, for example, if the signal levels are very low (for example, less that −20 dBFS); therefore, for the most accurate estimation, the process can continually monitor the incoming signal until its amplitude is above a predetermined threshold. The single amplitude and frequency content can be measured, in one embodiment, with a FFT analysis of the input signals. A 512 point FFT can be performed with a microprocessor or embedded processor, and the magnitude of each FFT bin can be measured to determine the signal amplitude at various frequencies.

If the signal amplitude and frequency content is suitable, then the nonlinear distortion levels can be estimated (step 715). Optionally, the current estimations can be compared (step 720) to the factory calibration values stored in the catalog to insure that the current estimations are not significantly different. A significant difference can indicate a component failure, a signal glitch, or other anomaly, whereby a determination can be made that the current estimation is not valid. If the current estimation is determined to be valid, then the new estimates are stored (step 725) in the catalog.

Based on the current operating conditions, such as current operating mode 635, temperature 645, and time 655, the catalog is searched and any relevant measurements are extracted (step 730), such as measurements near the current system temperature or having the same operating ode. If any pertinent measurements are missing from the catalog, they can be interpolated or extrapolated (step 735). For example, the catalog may contain estimates at a few different frequencies, but the values at the remaining frequencies can be interpolated or extrapolated using, in at least one embodiment of the invention, spline data fitting algorithms. Based on the interpolated/extrapolated data, the optimal nonlinear distortion parameters 155 and 165 can be calculated (step 740). The parameters are used to update (step 745) the nonlinear models 210A-N and the gain/phase shift filters 220A-N

Referring back to FIG. 5, the controller 570 operates in the background to guide the processing of data with the spectrum analyzer 530, channelizer 540, transfer function calculator 550, coefficient generator 560, and nonlinear distortion cataloger 610. The controller 570 may be implemented with a digital signal processor or an embedded processor, for example, a field programmable gate array. A program instantiated in the controller 570 can implement the process 700 for heuristically compensating nonlinear distortion, as shown in FIG. 7.

The controller 570 may also implement an iterative optimization algorithm that repeatedly analyzes the error signal via the spectrum analyzer 530 to adjust the transfer function calculator 550 and the coefficient generator 560 for optimal performance. Standard iterative optimization algorithms, such as a binary search or Newton's Method may be used. For example, the filter coefficients 165 calculated by the coefficient generator 560 may be iteratively optimized to adjust the amplitude and/or phase shift of the predistortion signal such that the selected distortion components are nulled in the output.

In general, the power specification for an amplifier is chosen based on the type of signal to be transmitted (e.g., the peak-to-average signal level), the desired power output, and the desired dynamic range (e.g., third-order intercept point (IP3), 1 dB compression point). Conventional amplifiers are often backed off such that the peak power does not exceed the 1 dB compression point (typically, amplifiers are backed off even a few dB more to insure signals remain in the linear operating region). For example, a typical communications signal may have a peak-to-average ratio of 9 dB, so the amplifier may be backed off by approximately 12 dB below its 1 dB compression point to insure linear amplification. The present invention can be used to increase the 1 dB compression point by 3 to 6 dB, which allows the back-off to be reduced commensurately. This corresponds to reducing the necessary power rating for the amplifier by one-half to one-quarter, which significantly improves the amplifier efficiency (i.e., as the back-off decreases, the efficiency increases). Moreover, the present invention provides a 35 to 40 dB improvement to the spurious free dynamic range (SFDR). Conventional linearization techniques only provide a 10 dB improvement.

Standard, commercially-available field programmable gate array (FPGA) chips are capable of digital signal processing at approximately 400 MHz data rate. For data rates higher than 400 MHz, the processing can be transformed into simultaneous parallel channels, for example, using polyphase filtering structures. For example, 4 GHz data can be demultiplexed into 16 lower data rate (250 MHz each) channels for real-time parallel processing in a standard FPGA (providing an instantaneous bandwidth of 2 GHz).

Linearity compensation techniques are described in commonly-owned U.S. Pat. No. 6,198,416, and calibration techniques are described in U.S. patent application Ser. No. 12/112,380, entitled “Adaptive Mismatch Compensators and Methods for Mismatch Compensation,” the disclosures of which are incorporated by reference herein in their entirety.

The following describes additional layers of processing that can be added to the linearization techniques described above.

FIG. 8 illustrates a block diagram of a linearity compensator 800 according to an embodiment of the invention. The linearity compensator 800 operates on an analog-to-digital converter (ADC) 810 and comprises the distortion compensator 100 and one or more factored Volterra compensators, which may include a second-order factored Volterra compensator 820, a third-order factored Volterra compensator 830, and additional higher-order factored Volterra compensators 840. Output of the ADC 810 is coupled to the inputs of the distortion compensator 100, the second-order factored Volterra compensator 820, the third-order factored Volterra compensator 830, and the additional higher-order factored Volterra compensators 840. The outputs of the distortion compensator 100, the second-order factored Volterra compensator 820, the third-order factored Volterra compensator 830, and the additional higher-order factored Volterra compensators 840 are coupled to an adder as shown. Volterra (or Volterra series) is a model for non-linear behavior where the output of the nonlinear system depends on the input to the system at theoretically all other times, thereby modeling intermodulation distortion of the system. Inclusion of the factored Volterra distortion compensators 820-840 improves linearization processing performance while significantly reducing the computational complexity compared to a traditional Volterra-based compensator without using the distortion compensator 100. In an embodiment of the invention, the factored Volterra distortion compensators 820-840 can be bypassed where only the distortion compensator 100 is implemented as described in FIGS. 1-7.

FIG. 9 illustrates the structure of the second-order factored Volterra compensator 820 according to an embodiment of the invention. This second-order factored Volterra compensator 820 comprises a number, N, of linear filters 910A-N, each of which has a corresponding squaring function 920A-N. The outputs of the squaring functions 920A-N are added together to produce the compensation signal output of the second-order factored Volterra compensator 820. In an embodiment of the invention, filter 910A is a 4-tap FIR filter. Accordingly, N is equal to 4, i.e., there are four filters 910A-D. The N linear filters 910A-N correspond to the dominant factors of the Volterra series; the process of measuring and factoring the Volterra series is described below.

FIG. 10 illustrates the second-order factored Volterra compensator operating in a parallel polyphase configuration 1000 according to an exemplary embodiment of the invention. Here, the ADC 810 is operating at 25 GHz. A 1:64 demultiplexer 1010 (“DEMUX”) is coupled to the output of the ADC 810 and a polyphase second-order factored Volterra compensator 1100. Accordingly, the polyphase second-order factored Volterra compensator 1100 operates on 64 demultiplexed ADC output signals, referred in the figure as C₀-C₆₃ at a relatively low data rate of 390.625 MHz. The demultiplexed compensation signal outputs from the polyphase second-order factored Volterra compensator 1100 are referred in the figure as F₀-F₆₃. Without this parallelization of the processing into numerous lower data rate paths, the extremely fast data rate would be beyond the capabilities of realizable hardware, such as digital signal processors (DSP), field programmable gate arrays (FPGA), or application specific integrated circuits (ASIC). Such parallelization may be implemented, for example, as polyphase finite impulse response (FIR) filters, the implementation of which is readily apparent to one of ordinary skill in the art. Likewise, a polyphase third-order factored Volterra compensator (not shown) operates on 64 demultiplexed ADC output signals, referred in the figure as C₀-C₆₃ at a data rate of 390.625 MHz.

FIG. 11 illustrates filters 910A-D implemented in a parallel polyphase configuration 1100 according to an embodiment of the invention. Each filter 910A-D is implemented using 64 parallel polyphase filters, each operating at 1/64 the system data rate. Here, because each filter 910A-D uses 4 taps (or multipliers) and there are 64 parallel polyphase filters, a total of 256 active multipliers are employed in the polyphase version of each filter 910A-D. This is repeated for each of the four filters 910A-D, thus using a total of 1,024 multipliers. The squaring functions 920A-N implemented in the parallel polyphase configuration use a total of 64 multipliers, one for each output F₀-F₆₃. Therefore, a total of 1,088 multipliers (1,024 plus 64) are needed for the parallel polyphase configuration 1100 of the second-order factored Volterra compensator 820

FIG. 12 illustrates the structure of the third-order factored Volterra compensator 830, according to an embodiment of the invention. The third-order factored Volterra compensator 830 comprises a number, N, of linear filters 1210A-N, the outputs of which are coupled to a 3^(rd) order polynomial function 1220. In an embodiment of the invention, filter 1210A is a 4-tap FIR filter. Accordingly, N is equal to 4, i.e., four filters 1210A-D. The 3^(rd) order polynomial function 1220 has 16 outputs which are coupled to an adder (not shown) to produce the compensation signal output of the third-order factored Volterra compensator 830. The filters 1210A-N can be implemented in a similar parallel configuration to 1100 in FIG. 11, where the squaring functions are replaced with 3^(rd) order polynomial functions. This structure of the third-order factored Volterra compensator 830 can be extended to higher order factored Volterra compensators by using higher order polynomial functions in place of the 3^(rd) order polynomial function 1220, the process of which is apparent to one of ordinary skill in the art.

Linear filters 910A-N and 1210A-N implement factored Volterra kernels which are measured via harmonic probing. Harmonic probing is implemented by injecting known multi-tone test signals into the ADC input. In a preferred embodiment of the invention, the number of tones needed is greater than or equal to the order of the kernel (e.g., greater than or equal to 3-tones for a third-order kernel). The frequency of each test tone is selected to minimize (or eliminate) overlapping components, e.g., harmonics do not have the same frequency as the fundamental tone. Overlapping measurements can be discarded so that only non-overlapping measurements are employed. More tones can be used than required for fewer test signals (e.g., using 3-tone test signals to measure second- and third-order kernels), with the same constraints on overlapping components. Test frequencies are selected to be centered on FFT bins (e.g., prime-coherent sampling) and all combinations are chosen to cover the full N-dimensional frequency space (e.g., for 3-tone combinations, each tone covers the desired frequency range and are calculated, for example, with embedded for-loops for each tone). Symmetry of the Volterra kernels reduces the required number of tone combinations. The amplitude of the test signals are preferably set to near full scale, i.e., saturation, of system where the distortion is most prominent.

In an alternative embodiment of the invention, test frequencies are selected randomly for uniform coverage over the N-dimensional frequency space. Test signals with too many overlapping components may be discarded. In order to prevent introduction of external non-linear distortion, high-linearity signal combiners are employed and the signal generator outputs' harmonic distortion is filtered via low-pass or band-pass filters.

A frequency domain Volterra kernel is analogous to a multi-dimensional linear filter where the input is

X _(n) [k ₁ ,k ₂ , . . . ,k _(n) ]=X[k ₂ ]*X[k ₂ ] . . . *X[k _(n)]

and the output is an n-dimensional convolution of input

X _(n)(z ₁ ,z ₂ , . . . ,z _(n))=X _(n)(z ₁)X _(n)(z ₂) . . . X _(n)(z _(n))

with n-dimensional Volterra kernel H_(n)(z₁, z₂, . . . , z_(n)). The one-dimensional output in the time domain is the diagonal of n-dimensional output y_(n)(k, k, . . . k).

Since multi-tone sinusoidal inputs in time-domain correspond to dirac delta functions in frequency-domain, the frequency response of the Volterra kernel H_(n)(z₁, z₂, . . . z_(n)) is effectively sampled at the multi-tone frequencies by using multi-tone sinusoidal test signals. The complex frequency response of the one-dimensional output is calculated, for example, with a one-dimensional FFT and evaluated at the frequencies of the harmonics and intermodulation distortion components (i.e., the “sum” frequencies). Care must be taken to appropriately scale the amplitude of the complex frequency response if some of the frequencies in the multi-tone sinusoidal test signal are repeated. This process is repeated for multiple multi-tone sinusoidal such that the Volterra kernel H_(n)(z₁, z₂, . . . z_(n)) is cataloged for many frequencies over the desired bandwidth.

In one alternate embodiment of the invention, the factored Volterra compensator is operated at a sample rate that is commensurate with the bandwidth of the nonlinear distortion that is being modeled. For example, if the maximum frequency in the desired band is f₁ and the maximum order of the compensator is 3, then the maximum frequency of the nonlinear distortion is the third harmonic at 3f₁. Therefore, the factored Volterra compensator can be operated at a minimum sample rate of 6f₁ to properly resolve the bandwidth of the modeled distortion components according to the Nyquist theorem. Bandpass interpolation can be used to appropriately increase the same rate of the factored Volterra compensator. Increasing the sample rate of the compensator has the effect of only utilizing particular subbands of the Volterra kernel since nonlinear distortion components will not fall at all possible frequencies. The filter design algorithm (described below) can be weighted to ignore these “do not care” bands and improve performance in the active subbands.

The efficiency of the measurement of the Volterra kernels is greatly simplified due to the symmetry of the Volterra kernels, h_(n)[k₁, k₂, . . . k_(n)]. The kernels h_(n)[k₁, k₂, . . . k_(n)] are equal for all permutations of k₁, k₂, . . . k_(n). This leads to very symmetric kernels in both the time-domain and the frequency-domain and super-symmetric factorization, which is described below.

To account for parameters that may drift over time and temperature, the device or system being compensated can be taken offline periodically for recalibration using the multi-tone harmonic probing approach discussed above.

In an alternate embodiment of the invention, the system can be adaptively calibrated in accordance with methods described above. In addition, correlation techniques can be used to estimate Volterra kernels for arbitrary inputs. Orthogonal factorizations of the Volterra kernels can be used for statistical independence, thereby simplifying the correlation measurements. The adaptive calibration can intelligently update factory calibration measurements and combine measurements over time, temperature, signal content, etc. as described above.

In a preferred embodiment of the invention, the calculation of Volterra kernel coefficients is performed with an over-constrained linear least mean squares filter design with filter weighting. A matrix of Volterra kernel coefficients is scaled by the measured complex frequency response. The real and complex parts are evaluated and subtracted from the measurement. This process is repeated for at least as many measurements as there are unique Volterra kernel coefficients. The kernel coefficients are calculated with an over-constrained least-squares solution with optional weighting of the form

WAx=Wb,

Where W is the weighting function (to optionally weight certain measurements or frequencies), x is a vector of estimates of the time-domain Volterra kernel, b is a vector of real and imaginary frequency response measurements, and A is a matrix corresponding to the frequencies of each measurement. The solution for the optimal Volterra kernel x via the over-constrained least-squares problem is

x=inv(A′W′WA)A′W′Wb)

This represents a very efficient, non-iterative solution via matrix algebra. Weighting is used to help ignore bands where signals will never fall. This approach can be extended to a MINIMAX (minimize maximum error signal) by iteratively adjusting the weighting of each measured frequency by a factor proportional to its amplitude. Therefore, larger amplitude error signals will get a higher weighting. The iteration is continued until it converges to the MINIMAX solution.

In an alternate embodiment of the invention, the calculation of Volterra kernel coefficients is performed via interpolation and extrapolation to uniform frequencies such that the computationally-efficient inverse Fourier transform can be used to estimate the kernel. Data at measured frequencies can be interpolated and extrapolate (e.g., using a cubic spline) to a set of uniformly-spaced frequencies corresponding to those of a Fast Fourier Transform (FFT). The inverse FFT translates these frequency-domain measurements to time-domain Volterra kernels. Forcing conjugate symmetry in the FFT measurements insures that real-valued (not imaginary) coefficients will be calculated for the Volterra kernels. Most stable, realistic systems exhibit fading memory where the amplitudes of the Volterra kernels decrease to negligible levels as the Volterra matrix gets larger. This inverse FFT filter design method can be used to easily calculate very large Volterra matrices, much longer than the anticipated fading memory size. A large inverse FFT also avoids frequency aliasing which could otherwise decrease the accuracy of the calculated Volterra kernel coefficients. A subset of the Volterra kernel can be selected by windowing the matrix centered on the maximum magnitude kernel coefficients. However, this method does not have a frequency weighting capability, so it may not be the optimal solution for a different error criterion. In that case, it may be used as starting point for other optimizations, such as an iterative solution.

In another alternate embodiment of the invention, the calculation of Volterra kernel coefficients is performed via an iterative linear or non-linear optimization. Volterra kernel coefficients can be iteratively adjusted and compared in a mean squares sense to the measured data, and the iteration continues until it converges to a solution.

For higher-order implementations of the factored Volterra compensator (i.e., order greater than 3), it becomes necessary to discriminate between overlapping kernels. Higher-order kernels have overlapping frequency components with lower-order kernels (e.g., some fifth-order intermodulation components are the same frequency as the third-order intermodulation components). In a preferred embodiment of the invention, the higher-order kernels can be measured at non-overlapping frequencies, followed by measuring the lower-order kernels by subtracting out overlapping components. In an alternate embodiment, kernels can be measured at multiple signal amplitudes and the different orders can be discriminated by their variations in amplitude (e.g., if the input amplitude is lowered by 1 dB, then the third-order components are reduced by 3 dB and the fifth-order components are reduced by 5 dB).

While it is possible to implement Volterra kernels as brute force multiplication and addition of all components, this is a very inefficient method requiring impractical amounts of processing resources, large size, high power consumption, and high cost. In a preferred embodiment of the invention, the Volterra kernel is factored into dominant components and a very accurate but efficient implementation of the kernel is performed by implementing only the most significant components. Insignificant factors can be ignored to significantly reduce the size with negligible impact on the performance. A factorization (also called decomposition) method called Singular Value Decomposition (SVD) is effective for second-order Volterra kernels, and a factorization method called Tucker Decompositions is effective for higher-order Volterra kernels (i.e., 3^(rd) order and greater). As discussed above, symmetry of the Volterra kernel leads to super-symmetric factorizations, which are extremely efficient. Note that these factorizations are exact implementations of the Volterra kernel if all factorization components are implemented. However, the benefit of the factorizations is the ability to decompose the Volterra kernels into its factors, rank the factors by the significance of their contribution, and elect to implement the most significant factors necessary for the desired level of performance.

In a preferred embodiment of the invention, the Volterra kernels are implemented as a parallel filter bank followed by a memoryless nonlinear combination of polynomials. If the memory of Volterra kernel is L (i.e., h_(n) is an L×L×L× . . . tensor), then there are at most L filters in the parallel filter bank, each with L coefficients. The maximum number of memoryless nonlinear combination is L^(k) (L=memory, k=Volterra order). The coefficients are quantized for implementation in efficient fixed-point arithmetic. In an alternate embodiment of the invention, the implementation can be floating point arithmetic for improved accuracy but requiring more processing resources.

In an alternate embodiment of the invention, the Volterra kernels can be factored by implementing only the dominant diagonals of Volterra kernel matrices. Often, the dominant energy in the Volterra kernel is concentrated on a few diagonals, so this factorization method can provide high accuracy with low computational complexity.

Once the Volterra kernels have been factored, they are rank ordered according to their significance (e.g., their singular values, Tucker factors, or other measurement of the relative energy in the Volterra kernel). Factored components are progressively included in the implementation until a desired level of performance or computational complexity limit has been reached.

Once the Volterra kernels have been decomposed into their dominant factors, the compensation system is implemented by negating the Volterra kernels above the first order. The first order term passes the fundamental signal through the compensator in phase, and the kernels above the first order are negated such that they are subtracted from the output, thereby canceling the nonlinear distortion.

FIG. 13 depicts an alternate embodiment of the invention, where a subset of the full Volterra kernel is implemented in a Linear-Nonlinear-Linear (LNL) Hammerstein-Wiener non-linear error compensation model. This model exercises much more of Volterra kernel than just Hammerstein or Wiener model by itself and is much more efficient than full Volterra implementation. The LNL compensator 1300 is an alternate embodiment of the distortion compensator 100 shown in FIG. 1. Referring back to FIG. 13, the LNL compensator comprises a first linear filter 105 which provides a fundamental signal path and a compensator signal path. The compensator signal path includes a distortion transfer function 130, which is comprised of a nonlinear polynomial 133 and a second linear filter 137. The fundamental signal path and the compensator signal path are summed via an adder to provide the compensator output, labeled z[n] in the figure. Optionally, the compensator path may include a bandpass interpolator, which is comprised of an upsampler 110 and a bandpass filter 137. The bandpass interpolator is used to resolve super-Nyquist distortion components. This option also includes a downsampler 140 to re-establish the original sample rate.

Measurement of the LNL kernel coefficients (i.e., the filter coefficients for linear filters 105 and 137) is essentially identical to that of the full Volterra kernel via harmonic probing via injection of known multi-tone test signals. Unlike the full Volterra harmonic probing, 1-tone and 2-tone measurements are sufficient for measuring higher order LNL kernels (i.e., measuring the LNL coefficients does not require at least as many tones as the order of kernel). The same requirements for choosing non-overlapping test frequencies apply. The same methods for measuring the complex frequency response at the sum frequencies using FFTs and cataloging the frequency response for many frequencies over the desired bandwidth applies. Similarly, in one alternate embodiment of the invention, the LNL compensator is operated at a sample rate that is commensurate with the bandwidth of the nonlinear distortion that is being modeled.

The equivalent Volterra kernel for the LNL compensator has the same Volterra symmetry, i.e., the equivalent kernels h_(n)[k₁, k₂, . . . k_(n)] are equal for all permutations of k₁, k₂, . . . k_(n). This leads to very symmetric kernels in both the time-domain and the frequency-domain and super-symmetric factorization. The equivalent Volterra kernel for the LNL compensator (i.e., the system with two FIR filters a[n] 105 and b[n] 137 separated by a polynomial operator, such as squaring or cubing) is the two-dimensional convolution of (a[n]′*a[n]) and a diagonal matrix that is comprised of the coefficients of b[n] as its diagonal components. The equivalent Volterra kernel frequency response is A(z₁)A(z₂)B(z₁z2) for second order, A(z₁)A(z₂)A(z₃)B(z₁z₂z₃) for third-order, and so on for higher order.

One method for measuring the coefficients of the first FIR filter a[n] 105 is with the injection of single-tone test signals with known frequency and phase. The transfer function (amplitude and phase) is measured by comparing the output amplitude and phase to the input amplitude and phase. A set of linear filter coefficients are fit to the measured transfer function using standard filter design techniques.

The same methods for recalibration of the Volterra kernel apply to the LNL compensator (i.e. periodic recalibration or adaptive background calibration). Several methods are available for the calculation of the filter coefficients for the LNL compensator. One method is an iterative non-linear optimization that uses a least mean squares curve fitting of the system to measured data (real and complex parts for both filters simultaneously). This method can be augmented with iterative filter weighting, for example, for MINIMAX optimizations. Another method optimizes the second filter, b[n], using previously described MINIMAX techniques and given a previous measurement of the first filter (a[n]). Another method is a mathematical projection of the full Volterra kernel into the LNL filter space. An unconstrained Volterra kernel is calculated with methods described above, and this solution is optimally projected given the available computational constraints of the LNL system.

As noted above, the LNL compensator is implemented as two FIR filters 105 and 107 with a memoryless polynomial nonlinearity 132 between the filters. The above filter design techniques are used to calculate the required filter coefficients, which can be truncated or windowed to extract the dominant coefficients. The coefficients can be quantized for an efficient fixed-point arithmetic implementation. The gain and linear delay of the system can be seamlessly shifted between the two filters with no practical effect on the overall response of the system.

The present invention is applicable to a wide range of military and commercial applications including, but not limited to: advanced radar systems; software-defined radios; multi-beam adaptive digital beamforming array transceivers, smart radios for wireless communications (terrestrial and satellite); wideband electronic warfare transceivers; general test equipment such as oscilloscopes, spectrum analyzers, and network analyzers; special test equipment, wide bandwidth modems, anti-jam global positioning system (GPS) receivers, and active radar for Earth science measurements.

The invention has been described herein using specific embodiments for the purposes of illustration only. It will be readily apparent to one of ordinary skill in the art, however, that the principles of the invention can be embodied in other ways. Therefore, the invention should not be regarded as being limited in scope to the specific embodiments disclosed herein, but instead as being fully commensurate in scope with the following claims. 

I claim:
 1. A linearity compensator for removing nonlinear distortion introduced by an electronic system, the linearity compensator comprising: a distortion compensator; and one or more factored Volterra compensators.
 2. The linearity compensator of claim 1, wherein the one or more factored Volterra compensators comprise a second-order factored Volterra compensator.
 3. The linearity compensator of claim 1, wherein the one or more factored Volterra compensators comprise a third-order factored Volterra compensator.
 4. The linearity compensator of claim 3, further comprising additional higher-order factored Volterra compensators, wherein the additional higher-order factored Volterra compensators are fourth or higher-order factored Volterra compensators.
 5. The linearity compensator of claim 2, wherein the second-order factored Volterra compensator comprises a number of parallel finite impulse response filters and a number of memoryless second order polynomials.
 6. The linearity compensator of claim 5, wherein the second-order factored Volterra compensator is implemented in a parallel polyphase configuration.
 7. The linearity compensator of claim 3, wherein the third-order factored Volterra compensator comprises a number of parallel finite impulse response filters and a number of memoryless third order polynomials.
 8. The linearity compensator of claim 7, wherein the third-order factored Volterra compensator is implemented in a parallel polyphase configuration.
 9. The linearity compensator of claim 1, wherein the electronic system comprises an analog-to-digital converter or a power amplifier.
 10. The linearity compensator of claim 1, wherein the one or more factored Volterra compensators comprises a number of linear filters implementing factored Volterra kernels.
 11. The linearity compensator of claim 10, wherein the factored Volterra kernels are measured via harmonic probing comprising multi-tone test signals injected into the input of the electronic system.
 12. The linearity compensator of claim 11, wherein a number tones in the multi-tone test signals is greater than or equal to the order of the corresponding factored Volterra kernel.
 13. The linearity compensator of claim 10, wherein the factored Volterra kernels are measured via harmonic probing comprising random test tone signals with non-overlapping distortion components.
 14. The linearity compensator of claim 1, wherein the linearity compensator comprises an adaptive nonlinear distortion estimator.
 15. The linearity compensator of claim 10, wherein the factored Volterra kernels are calculated using an over-constrained linear least mean square filter design.
 16. The linearity compensator of claim 10, wherein the factored Volterra kernels are calculated using interpolation and extrapolation and inverse fast Fourier transforms.
 17. The linearity compensator of claim 1, wherein the distortion compensator comprises: a first linear filter providing a fundamental signal path and a compensator signal path; a distortion transfer function operating on the compensator signal path, wherein the distortion transfer function comprises a nonlinear polynomial and a second linear filter; and an adder for summing the fundamental signal path and the compensator signal path.
 18. The linearity compensator of claim 17, further comprising a bandpass interpolator operating on the compensator signal path.
 19. The linearity compensator of claim 18, wherein the bandpass interpolator comprises an upsampler and a bandpass filter.
 20. The linearity compensator of claim 18, further comprising a downsampler operating on the compensator signal path.
 21. The linearity compensator of claim 17, wherein filter coefficients of the first linear filter are measured by fitting the filter coefficients to a coherently measured amplitude and phase of a transfer function of the distortion transfer function.
 22. The linearity compensator of claim 17, wherein the distortion transfer function further comprises an adaptive nonlinear distortion estimator.
 23. The linearity compensator of claim 17, wherein filter coefficients of the first linear filter are iteratively optimized.
 24. The linearity compensator of claim 17, wherein filter coefficients of the first linear filter are optimized using MINIMAX optimization.
 25. The linearity compensator of claim 17, wherein filter coefficients of the second linear filter are optimized using MINIMAX optimization
 26. The linearity compensator of claim 1, wherein the factored Volterra compensators comprise a subset of all the factors. 